Electronic Structure of Lanthanum Hydrides 
with Switchable Optical Properties 
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Recent dramatic changes in the optical properties of LaHb+i and YH2+1 films discovered by Huib- 
erts et al. suggest their electronic structure is described best by a local model. Electron correlation 
is important in i?~-centers and in explaining the transparent insulating behavior of LaH3. The 
metal-insulator transition at x ~ 0.8 takes place in a band of highly localized states centered on the 
JLvacancies in the LaH3 structure. 



Recently Huiberts et. al. pj reported dramatic changes 
in the optical properties of lanthanum and yttrium hy- 
dride films with changing hydrogen content, e.g. films 
that switch from a shiny mirror to a yellow transpar- 
ent window in a time of order seconds. Although many 
metal-insulator transitions are known, it is very un- 
usual and potentially of technological importance, that 
the transition leads to spectacular effects in the visible. 
Huiberts et. al. 0] point out that the electronic structure 
that underlies this behavior is poorly understood - in- 
deed standard LDA-calculations do not predict a metal- 
insulator transition at all [^],[| . In this letter we examine 
the role of electron correlation in these hydrides and ar- 
gue it justifies a local description that derives from an 
ionic starting point. 

The structural changes are especially small in the 
LaH2+x films. The La-atoms always form a fee lattice 
with two Lf-atoms always occupying tetrahedrally coor- 
dinated sites. As x changes from to I, the octahedrally 
coordinated site goes from empty to fully occupied and 
a good metal (LaLf 2 ) evolves to a transparent insula- 
tor (LaHs). The metal- insulator transition in the d.c. 
conductivity JO occurs at an intermediate concentration, 
x ~ 0.8, far removed from a lightly doped semiconduc- 
tor. This in turn suggests highly localized states asso- 
ciated with H- vacancies in LaH^. The key challenge is 
to understand why the structural changes with changing 
hydrogen content are remarkably small, but the changes 
in spectral properties are dramatic at energy scales up to 
a visible and secondly the origin of the highly localized 
states at intermediate concentrations. We shall address 
these issues in turn. 

We have performed self-consistent LDA calculations on 
LaH2 and LaH^. In agreement with previous results |^] 
there are two sets of states - low energy primarily ls- 
Lf-states and higher energy states of mainly 5d-La char- 
acter. These results suggest assigning a formal valence 
H~. In LaH2 the Lf~-bands are filled leaving lel/f.u. 
in the 5d-La conduction bands leading to metallic be- 
haviour 0]. In LaHz the Lf~-bands can hold all 6eZ/ f.u. 



but in the LDA calculations there is overlap between the 
H ~ and bd-La bands leading to a semimetal rather than 
the observed transparent insulator. 

The LDA calculations give an interesting insight why 
atomic H is easily incorporated in this structure. The 
net charge change between LaHi and LaH^ is found to 
be remarkably small. The charge at the octahedral site 
(within a sphere of 2.533 atomic unit (a.u.) radius) in 
LaHz is 1.588 el., while the charge in an empty sphere of 
the same radius in LaH 2 is 0.568 el., so to a very good ac- 
curacy, LaHi+x can be interpreted as neutral Lf-atoms 
moving into or out of the octahedral interstitial sites, 
with negligible changes in the lattice structure. This 
charge distribution of LaH^ is not very different in an 
ionic picture due to a large radius (~ 2a. u.) of the outer 
electron in the bound state of the Lf~-ion. 

The Lf~-ion is a difficult case for LDA and a careful 
treatment of the correlation between the two electrons is 
required to obtain the binding energy (~ 0.7 eV). 
This led us to examine the effect of correlations on the 
width of the Lf~-bands. We start by considering the 
hopping integral t(d) for an electron between H- atoms 
with separation d (in a.u.), ignoring for the moment the 
La-ions. This can be obtained at once from the lowest 
energy levels of the iLJ-ion with inter-proton distance, 
d. H2 is a bound negative ion for d > 3 - the range of 
interest here. 

Let ?/>j(l,2) be the groundstate wavefunction of two 
electrons in H~, and <f>i(l) be the groundstate of hy- 
drogen atom. The three electron states of the 
with odd(-) and even(+) parities may be constructed 
using the single site (i and j) states: \P± = ± 
where $j = A[ipi(l,2)<f)j(3)], and A is an anti- 
symmetrizer. The corresponding energies are given by 
E± = ( v E'±|/i| v E , ±)/( v I / ±| v I / ±), where h is the Hamiltonian 
for the iL7 -system, including the kinetic and Coulomb 
energies M. We choose a variational state of Chan- 
drasekhar [pj for H~: 
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where r are the electron radial coordinates, \ is the spinor 
for a singlet, and a = 1.075, b = 0.478, c = 0.312 in a.u.. 
This simple wavefunction describes two intraatomic cor- 
related electrons well, and its energy is very close to the 
best estimate of Hylleraas || involving 24 variational pa- 
rameters. The states for if^ thus constructed become 
exact in the limit d > 1. The effective hopping inte- 
gral for an electron between two -ff-sites can be obtained, 
t(d) = (E-—E + )/2. At d = 4, the groundstate energy es- 
timated in this method is E- = —1.032 a.u., very close to 
the best available result of -1.034 a.u. for using more 
complicated variational method B. This suggests that 
the groundstate of at these separations is very well 
described by a linear combination of H~ and if-states, 
and justifies our approach. 

The two shortest inter-hydrogen distances in LaH^ are 
those between n.n. tetra- and octa-sites (denoted as H tc t 
and -ffoct) and n.n. tetra-sites which denote as t 2 and 
t\ respectively. These distances are ^/3ao/4 and ao/2 
(ao: the lattice constant) respectively (4.58 and 5.29 a.u.) 
giving values t 2 = - 0.74 eV and t\ = -0.54 eV. The H~- 
bands in tight binding representation are the eigenstates 
of the Hamiltonian 

H th = £V(£)M(£) C (£), 

k 

where c\k) is a 3-component vector, representing the 3 
ff-states in a unit cell, and 

/ e t ha(k) t 2 f3(k) \ 
M(k) = t ia (k) e t t 2 /3*(k) . 
\t 2 (3*(k) t 2 (3{k) e ) 

with 

a (£) = X! 2cos(k T /2), 

T—x^y^z 

(Tl ,(T2=±1 

Note we now use the lattice parameter ao(5.60A) as the 
length unit. The atomic energy levels at H tet and H oct 
(et and £o) we estimate using the LDA calculations as 
follows. At the high symmetry T point, one can identify 
in the LDA a 4 by 4 matrix corresponding to the three H 
and one La-6s states in a unit cell of the fee lattice. We 
can map this matrix onto a 3 by 3 submatrix describing 
the effective three 7J-bands with the same three lower 
eigenenergies. The highest energy band, representing the 
La-6s state, is pushed well above the bottom of the La — 
5d conduction band due to the strong hybridization, and 
can be projected out. Incorporating the binding state 
energy of a single H~-ion, which is 0.7 eV, we obtain 
e t = —3.9 eV, and £o = —3.3 eV. 



In Fig. la the resulting energy bands are plotted. We 
find a small gap, which is not sufficient for a transparent 
insulator. However we have neglected the La 3+ -ions in 
estimating (ti, t 2 ). In an ionic picture the La 3+ -ions gen- 
erate a crystal field at the iJ-sites whose effect can be esti- 
mated as follows. Let the interaction between an electron 
in H2 and the other ions be h! = — J^fi Z(R)e 2 /\r — R\, 
with Z = 3 for La 3+ , and Z=-l for iJ~-ions, and R the 
ion position. The energies for of two H te t are mod- 
ified, due to the crystal field, to E'±, which is obtained 
by adding h! to h. The reduction of t\ due to the crystal 
field is found to be Sti = 0.20 eV, where the contribu- 
tion from within a cubic unit cell is 0.08 eV. The esti- 
mate for t 2 is more complicated because of the Madelung 
constants. Assuming the same percentage reduction, we 
estimate the hoppings in the presence of the crystal field, 
ti = -0.34 eV, i 2 = -0.47 eV. The estimate is based on an 
ideal ionic approach, and neglects the screening effect, so 
that the actual reduction of t's are surely smaller. With 
the new values (t\, t 2 ), we obtain strongly narrowed H~- 
bands (see Fig. 1) and a transparent insulating ground 
state. This idealized ionic approach overestimates the 
narrowing and neglects the effects of hybridization on 
the hopping integrals. These however are weak at the 
top of the _ff~-bands and are mainly important only at 
energies well below the Fermi energy. 

Turning to the metal-insulator transition in LaH 2+x as 
x increases, we can view it starting either at the metallic 
(x = 0) or insulating (x — 1) endpoints. In the former 
case, the introduction of a neutral 7?-atom into a H oct - 
site creates a s=l/2 magnetic impurity, which couples 
to the conduction 5d — La electron spins. The effective 
Hamiltonian is a Kondo model, 

n eS = J2 £ $) d l^ d k,a + J /~2^ g ^ 

k,cr • 

where i runs over all the occupied H oct , and S and s are 
the electron spins of the neutral H and of the conduction 
electron states. The latter represent linear superposition 
of the La — bd — e g states around three degenerate X\ 
points, which is of s — wave symmetry and can couple to 
the ls — H state. The antiferromagnetic coupling J ~ 0.7 
eV in free space, and may be renormalized by a numerical 
factor in hydrides. In this Kondo description conduction 
electrons are captured by the neutral iJ-atoms at -ff ct- 
sites to form tightly bound singlets and LaH^ is viewed 
as a Kondo insulator with a large band gap. 

Starting from insulating LaH^, removing neutral H- 
atoms causes vacancies at octa-sites (H^ ct ) which donate 
an electron to the conduction band. In a conventional 
semiconductor such as in the doped Si, the impurity 
state is described by an effective mass theory, and the 
result is a hydrogen-like bound state with a large effec- 
tive Bohr radius a* B of order of hundred A due to the 
light mass and the large dielectric constant. The critical 
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impurity concentration x c at which the system becomes 
metallic is given by Mott criterion x c ~ (l/a^) 3 , where 
I is the interatomic distance. Because of a* B 3> I, x c is 
very small (0.1 percent for Si). The vacancy state in 
LaH2+x is very different, however. Experimentally, it is 
found that the semiconducting states extend to x = 0.25 
for lanthanum hydrides, and to even larger value for yt- 
trium hydrides Below we shall show that the va- 
cancy state in lanthanum hydrides is an octahedral s-like 
La — e g state with extremely small size. 

The La-5d-e g and H oct electrons hybridize strongly to 
form bonding (mainly H) and anitbonding (mainly e g ) 
states. A vacancy of H oct breaks the bonds locally and 
the e g electron becomes locally non-bonding, which has 
much lower energy than the antibonding e g states away 
from the vacancy. Therefore, in addition to serving as 
a positive charge center as in conventional n-type semi- 
conductors, H^ct creates a potential well for the e g state 
electron. The latter is non-perturbative, and is responsi- 
ble to the unusual concentration-dependence of the semi- 
conductor. 

Consider a s— symmetric octahedral e g state (S-state 
hereafter) around a H^ ct as shown in Fig. 2. The H^ ct 
vacancy breaks the bonds on the octahedron and re- 
duces the anitbonding energy of the S'-state. This S'- 
state is very localized because the neighboring octahe- 
dral s— state is antibonding and has a much higher en- 
ergy. We can estimate the depth of the potential well, 
Vb, within a perturbation theory upto the second order 
in V s d(j, the hopping integral between an atomic ds z 2_ r 2 
and its neighboring H oct with the orbits towards each 
other. Let A be the atomic energy difference between 
La-5d and H oct , then V = -6V? da /A. This is 2/3 of the 
antibonding energy 9V^ da /A at the Xi point in LaH 3 . 
The latter is estimated to be ~ 6 eV from LDA po| , so 
that Vb ~ -4 eV. We may also compare this S state to 
the non-bonding t2 g state. The latter is the lowest en- 
ergy conduction state of LaH^, and would be a starting 
point for an effective mass theory in conventional semi- 
conductors. The S state has much better kinetic energy 
as evidenced from LDA for LaH2 [^,|io|,0. A more de- 
tailed calculation based on a tight-binding model of n.n. 
hopping between La sites estimates the kinetic energy 
gain of the S state over the t2 g state is ~ 0.9 eV. This 
is about to compensate the loss of the antibonding en- 
ergy of the S state with the neighboring occupied H oct 
atoms, which is 1.5V^ da / A ~ 1 eV. The S state has bet- 
ter Coulomb attraction to the vacancy, because its orbit 
is towards the central effective charge. Therefore the S 
state should have lower energy than the t2 g state, and is 
the lowest energy state for H oct vacancy. We emphasize 
that the potential well generated by the H oct vacancy 
is short-ranged, different from the long-range Coulomb 
force. The effective Bohr radius of the impurity state is 
only half of the lattice constant. 

The localized nature of the impurity state in LaH2+ x 



is consistent with the temperature dependence of the re- 
sistivity data. Shinar et. al. Q reported a phase tran- 
sition at T = 250K, above which the resistivity shows 
variable range hopping, below that temperature, it is co- 
herent state of the defect band. Localized states are a 
prerequisite for variable range hopping. The coherent 
state accompanies the onset of superlattice order in the 
HY ct vacancies. 

The precise nature of the transition in optical proper- 
ties is not well understood yet. Since the d.c. conductiv- 
ity is insulating only in the disordered state of H^ ct , and 
is metallic in the ordered one, we argue that the tran- 
sition should be a disorder driven Anderson transition, 
instead of a Mott type correlation driven transition. 

In conclusion, we have studied the electronic structure 
of LaH2+x starting from H~ -ionic picture emphasizing 
the electron correlation. Our calculations give an insula- 
tor for LaH%, in agrement with the recent experiments. 
We propose that a iJ-vacancy in LaH^ introduces a s- 
wave state centred at the vacant octa-site,with a highly 
localized form. This state of the vacancy explains why 
the semiconductor LaH2+ x to be so stable in a wide range 
of concentration x (0.8 < x < 1), and is consistent with 
the experimental observed variable range hopping mech- 
anism for the resisitivity in these materials. The recent 
experimental activity of Huiberts et. al. jjj has raised 
many interesting questions in study of metal hydrides. 
The precise nature and the criterion of the optical tran- 
sition from transparent to reflecting needs be examined. 
The role of quantum diffusion of the hydrogen at low 
temperatures, the order-disorder transitions, and the in- 
teractions between octahedral H- vacancies, are all inter- 
esting open questions. 
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FIG. 2. Diagrammatic illustration of the proposed impu- 
rity state in LaH^. The centre circle represents a missing 
H, forming an n-type impurity center. The surrounding or- 
bits represent phases for the local s-like octahedral La-5d-e s 
state. 
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FIG. 1. The valence bands of LaHz calculated from the 
single hole hopping in H~ . The dashed lines indicate the 
bottom of the conduction bands, a): without crystal fields; 
b): with the crystal fields 
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